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ABSTRACT 

The relic abundance of thermal dark matter particles is generally assumed to be inversely 
proportional to their annihilation rate, which is therefore constrained by the present matter 



density, {cr ann v) ~ 10 lb Vt 



cm 3 



sec . Here we point out that much lower values of {o~ a 



are possible for heavy dark matter candidates (mx ^ 10 TeV) that couple to other particle species 
through the electroweak force. With heavy dark matter particles present the early universe may 
evolve according to the following scenario. After an early entry into matter-dominated phase, dark 
matter particles form self-gravitating microhalos. Collisional interaction between dark matter 
particles and the surrounding radiation field eventually leads to microhalos gravothermal collapse 
and annihilation of most dark matter particles. For sufficiently heavy dark matter candidates 
(mx > 10 TeV) the universe can return to radiation-dominated phase before the nucleosynthesis 
and thereafter follow the "standard" scenario. 



Subject headings: cosmology: theory - early universe 

1. Introduction 

While there is a strong evidence for the ex- 
istence of non-baryonic dark matter, little is 
known about its origins. In the currently popu- 
lar scenario, dark matter particles and antipar- 
ticles start in thermal equilibrium with radia- 
tion. When temperature drops below the par- 
ticle mass, mx, their comoving density begins an 
exponential fall until the annihilation time-scale, 
tann = (n-x (o 'ann v )) > becomes larger than the 
age of the Universe. From this point the comoving 
density of dark matter particles is assumed frozen. 
For the relic abundance of dark matter particles to 
be consistent with the present value of 51dm their 
annihilation rate at the time of freeze-out must 
equal (Kolb & Turner 1989) 



(o-annV) ~ 3 • 10 



" 26 cm 3 



(1) 



Unless dark matter is of non-thermal origin 
(Chung, Kolb & Riotto 1998), lower values of 
(o~annV) would overdose the universe and there- 
fore can be excluded. 

If dark matter particles interact through any 
known force, during the decoupling epoch, when 



T ~ mx, their annihilation cross-section can not 
exceed a ann < a/T 2 ~ ajm\. Using this fact, 
Griest & Kamionkowski (1990) obtained from (Q} 
an upper limit on the mass of dark matter particles 



m x < 300 TeV. 



(2) 



However, in this paper we point out that dark 
matter candidates with higher mass and lower an- 
nihilation cross-section are not necessarily incon- 
sistent with the present value of fidm- Unlike the 
standard model, which assumes the dark matter 
comoving density freezes out soon after the tem- 
perature of the Universe drops below mx , we show 
that it is possible for it to drop again later, after 
dark matter particles form virialized halos. 

2. Dark matter annihilation inside halos 

Unlike the standard model, where the matter 
becomes dominant only around z ~ 10 4 , in a Uni- 
verse with low (cr ann v) and high mx the initial 
transition to matter domination happens much 
earlier. Consequently the growth of density per- 
turbations and the formation of the first gravita- 
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tionally bound halos would be likewise shifted to 
a much earlier epoch. 

The collapse of dark matter particles into halos 
gives a strong boost to the dark matter annihi- 
lation rate. Whereas previously the time-scale for 
particle annihilation was increasing faster than the 
age of the Universe, freeze-out of the local density 
inside gravitationally bound halos causes the re- 
versal of this trend. Once the Hubble time begins 
to catch up with the annihilation time-scale, the 
dark matter comoving density is no longer con- 
stant. 

Subsequent evolution of the Universe depends 
on the nature of the dark matter particles. For 
purely collisionless dark matter it can be shown 
(see Appendix A) that the comoving matter den- 
sity eventually approaches an asymptote, p rn (l + 
z)~ 3 oc t -1 / 4 . This moderate rate of decline turns 
out to be insufficient for the Universe to return to 
radiation dominated phase before the nucleosyn- 
thesis. However, annihilation inside halos can be 
greatly sped up if dark matter particles interact 
with the surrounding radiation (made up mostly 
of relativistic quarks and leptons). 

At present weakly interacting massive parti- 
cles (WIMPs) are the most popular candidates 
for the role of dark matter. The collisional cou- 
pling between WIMPs and other species is usu- 
ally assumed to be negligible. At our epoch this 
assumption is generally correct, as the Universe 
temperature is far below the mass scale of W ± 
and Z bosons which mediate the weak force. How- 
ever, in the early universe, when T > 300 GeV, 
the coupling becomes very efficient. At high tem- 
peratures WIMPs scattering cross-section rises to 
o elastic ~ a2 /T 2 , making the ratio between the 
relaxation time and the age of the Universe 

Trel m X + T 

t Hubble a z m P i 

where mpi is the Planck mass. For mx < 10 13 
TeV and 0.1 < T < 10 13 TeV this ratio drops 
below unity, implying that WIMPs kinetic energy 
is tightly coupled to the radiation. 

Collisional coupling between dark matter and 
radiation has a two-fold effect on the formation 
of self-gravitating dark matter halos, which is ex- 
pected to begin after the matter-radiation equal- 
ity. On small scales, T V i r i a i < T 7 , the energy 
transfer from radiation to dark matter prevents 



the collapse of the low-mass halos. On large scales, 
T V i r iai > T.y, the effect is the opposite. Since 
self-gravitating objects have negative heat capac- 
ity, the energy loss to the surrounding radiation 
field would lead dark matter particles to collapse 
deeper into the gravitational well, thereby raising 
the virial temperature. Thus the continuous en- 
ergy loss to the background radiation results in a 
gravothermal collapse. 

The time-scale for particle diffusion in self- 
gravitating objects, which in our case also sets the 
time-scale for gravothermal collapse of larger ha- 
los, is 



where a e i ~ a 2 /T 2 is the scattering cross-section 
for the weakly interacting particles. If, as gener- 
ally predicted by the inflation models, the density 
perturbation spectrum is scale-invariant, then the 
velocity dispersion of dark matter particles in ha- 
los with a size above the damping scales should 
be of order v ~ 0.1c and the virial temperature 
T ~ O.Olmx- Thus in large halos the gravother- 
mal collapse proceeds on a timescale 




If halos consist equally of dark matter particles 
and anti-particles (i.e. there is no dark matter 
asymmetry) at some point during the gravother- 
mal collapse nearly all of them may annihilate into 
radiation. Provided the mass of the particle is 
greater than ~ 10 TeV, halos collapse and the 
resulting conversion of dark matter to radiation 
would end before nucleosynthesis and thereafter 
the evolution of the Universe can proceed accord- 
ing to the "standard" scenario. 

3. Discussion 

It follows from our scenario that the existence of 
stable particles with high mass and low annihila- 
tion cross-section does not violate the constraints 
on the present matter density. If several species of 



During the epoch of reionization, 6 < z < 20, inverse 
Compton scatterings between electrons and the CMB pho- 
tons have a similar effect on baryonic halos. However, there 
the gravothermal collapse is generally avoided either by en- 
ergy injection from stars and miniquasars or by hydrogen 
recombination, which depletes the free electron abundance. 
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such particles existed in the early universe, each 
could have been eliminated in separate periods of 
matter domination followed by return to radiation 
domination. 

To satisfy the constraints on abundance of the 
light elements the last reversal to radiation dom- 
ination must have happened prior to the nucle- 
osynthesis. Nevertheless it is possible for the tran- 
sitional matter dominated phase to leave its im- 
print on the presently observed universe. Besides 
leaving a small remnant of heavy particles, which 
at present compose the dark matter, the transi- 
tional matter dominated phase might affect later 
universe by creating a population of primordial 
black holes and baryon asymmetry. 

The violation of CP-symmetry in the early uni- 
verse (Sakharov 1967) makes it natural to expect 
a small differential between the coupling strengths 
of particles and anti-particles. Inside collapsing 
halos the drag force produced on dark matter 
particles through elastic scatterings with the es- 
caping radiation may be different from the force 
on dark matter anti-particles. As a result, spa- 
tial segregation of dark matter particles and anti- 
particles would be produced. Similarly, a dif- 
ferential drag by the infalling dark matter par- 
ticles on the outflowing quarks and anti-quarks 
(or lcptons and anti-leptons) would produce lo- 
cal over/undcrabundances of baryons and lcptons. 
During halos gravothcrmal collapse the existence 
of such segregation would make it impossible for 
the entire halo to annihilate into radiation. There- 
fore the collapse of particles (or anti-particles) 
dominating at the halo center should lead to the 
formation of a black hole. 

Since black holes do not conserve baryon num- 
bers (Hawking 1974; Zeldovich 1976; Carr 1976), 
formation of a central black hole would convert a 
local particle asymmetry into a global one. De- 
pending on their initial mass, these primordial 
black holes may have evaporated long before the 
present epoch, possibly further amplifying the 
matter-antimatter asymmetry in the process (Dol- 
gov 1982). Alternatively, if their mass was suffi- 
ciently large, significant population of primordial 
black holes may still be found in our Universe. 



The author thanks E.W. Kolb, M.S. Turner and 



M. Beltran for stimulating discussions. 

Appendix A: Annihilation of collisionless 
particles in virialized halos 

We consider the evolution of matter den- 
sity field, assuming the annihilation cross-section 
scales as (a ann v) oc v~ x at non-relativistic veloci- 
ties. While it is common to assume x = 0, many 
new particle candidates with x = 1 have been 
proposed (e.g. Hisano et al. 2004; Profumo 2005; 
Lattanzi & Silk 2008). Here we shall assume that 
x can take any value in the physically plausible 
range < x < 2. 

After decoupling from radiation matter den- 
sity drops to a point when the annihilation time- 
scale, t ann — (n x (a 'annv))" 1 , becomes longer than 
the Hubble time, so that the comoving density is 
frozen. As long as the density field remains nearly 
homogeneous with the density and velocity dis- 
persion falling respectively as a~ 3 and a" 2 with 
the increasing scale factor, a, t ann /tuubbie contin- 
ues to rise. However, this trend is reversed af- 
ter matter collapses into halos. Inside an isolated 
gravitationally bound halos both the local density 
and velocity dispersion remain roughly constant, 
freezing the growth of t ann until it becomes com- 
parable with t Hubble- At this point a halo would 
start losing mass at a significant rate via parti- 
cle annihilation. The subsequent halo evolution of 
can be described as following. 

The mass loss from particle annihilation and 
the resulting decrease of binding energy would 
cause the halo to expand, so that the annihila- 
tion time-scale tracks the age of the halo. Because 
the dynamic time is much shorter than t Hubble! 
the halo always remains close to dynamical equi- 
librium. Since the gravitational binding energy is 
proportional to E g oc M 2 /R, where M and R are 
the total mass and the radius of the halo, changing 
the mass and radius by dM and dR respectively, 
changes E g by 

dE ° = E » [ 2 W - -r) • (6) 

The change of halo kinetic energy, Ek , in the an- 
nihilation process depends on the velocity depen- 
dence of Tann and the particle velocity distribu- 
tion. Assuming Maxwell-Boltzmann distribution, 
we find that the average kinetic energy of an an- 
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nihilating particle is (1 — x/2)(mx (w 2 ) /2), where 
(v 2 ) is the local velocity dispersion. Thus 



dE k 
E k 



x\ dM 



A I 



(7) 



Combining the virial theorem, E g = —2Ek, with 
equations ([5]) and ([7]), we find that the radius 
of the halo increases with the falling mass as 
R oc M- 1+x / 2 . Further, by using the scalings 
p oc M/R 3 , v oc y/M JR and (pv~ x ) oc t~ x , we 
find that the total mass of the halo falls with time 
as M oc p, = -4/(16 + 2x - x 2 ). For the entire 
range < x < 2, p is confined to a narrow range 
-0.25 < x < -0.235. 

Modeling the Universe as an ensemble of iso- 
lated halos, we can follow the global evolution of 
matter and radiation density fields 
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d 2 a 
adt 2 

dp m 

dt 

dp-y 

~dt 



(p-m + 2p-y 



da p 
3 ^dt ~ J ' 



da 
adt 



p 
J* 



(8) 
(9) 
(10) 



where p 1 and p m are respectively, the mean densi- 
ties of radiation and matter and a is the expansion 
factor. Solving the equations, we find the asymp- 
totic solutions 



a oc t^l*, 

2 + bp + 2p 2 



Pm 
Pi 

For p R 



12nGt 2 
2{l + 2p) 

3p 

-1/4 we get fi m 



Pi 



-(2 + ^ 
SnGt 2 



(11) 
,(12) 

(13) 



4/7 and fi 7 « 3/7. 
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